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The quantized lateral motional states and the spin states of electrons trapped on the surface of 
superfluid helium have been proposed as basic building blocks of a scalable quantum computer. 

Circuit quantum electrodynamics (cQED) allows strong dipole coupling between electrons and a 
high-Q superconducting microwave resonator, enabling such sensitive detection and manipulation 
of electron degrees of freedom. Here we present the first realization of a hybrid circuit in which a 
large number of electrons are trapped on the surface of superfluid helium inside a coplanar waveguide 
resonator. The high finesse of the resonator allows us to observe large dispersive shifts that are many 
times the linewidth and make fast and sensitive measurements on the collective vibrational modes 
of the electron ensemble, as well as the superfluid helium film underneath. Furthermore, a large 
ensemble coupling is observed in the dispersive regime during experiment, and it shows excellent 
agreement with our numeric model. The coupling strength of the ensemble to the cavity is found 
to be > 1 MHz per electron, indicating the feasibility of achieving single electron strong coupling. 


INTRODUCTION 

Electrons on helium are a promising resource for quan¬ 
tum optics and quantum computing[ - ]. They form an 
extremely clean two dimensional electron gas [5], as evi¬ 
denced by a mobility exceeding 10 7 cm 2 /Vs [6, ']. and 
the electron spin coherence time is predicted to exceed 
10 3 s [1]. Electrons on helium have been used to study 
Wigner crystallization and quantum melting [8-10]. Re¬ 
cent experiments employ them as a powerful probe to 
study the topological domain structures on the surface 
of superfluid helium 3 [11—13]. In addition, it is now 
possible to build mesoscopic structures such as quantum 
dots with one or a few electrons on helium, and single 
electron scale charge coupled devices [14-17]. However, 
performing quantum experiments in this fascinating sys¬ 
tem has lagged behind that in semiconducting 2D elec¬ 
tron gasses, such as GaAs, as traditional measurement 
techniques cannot be applied to electrons on helium. In 
particular, it is not possible to make direct Ohmic contact 
to the electron gas. Additionally, the largely unscreened 
electron-electron Coulomb force and a hydrostatic insta¬ 
bility of the system [18, 19] suppress the exchange inter¬ 
actions typically used in semiconductor spin qubits [20]. 

The circuit QED architecture [21, 22] offers a path 
to new experiments in the quantum regime as well 
as improving the sensitivity and bandwidth of existing 
measurements. In this hybrid approach, electrons are 
trapped above an on-chip superconducting microwave 
resonator. The presence of the electrons changes the ef¬ 
fective capacitance of the cavity, resulting in a dispersive 
shift of the cavity resonance frequency. In the strong dis¬ 
persive regime, the cavity frequency shift is larger than 
the cavity linewidth, and every photon measures the state 
of the electrons. Because the energy of a single photon in 
the cavity is higher than the thermal bath (fio; > k^T), it 


is possible to conduct quantum optics experiments at the 
single photon level. This dispersive measurement is con¬ 
ceptually similar to the Sommer-Tanner technique [23], 
but the use of resonant superconducting circuits at mi¬ 
crowave frequencies enables better impedance matching, 
resulting in faster and more sensitive measurements of 
small ensembles. Finally, the hybrid architecture allows 
one to leverage the substantial progress in superconduct¬ 
ing circuits over the past decade [24, 25]. 

In this letter, we report the first implementation of a 
circuit QED architecture with electrons on helium. We 
show lithographic control and sub-nanometer measure¬ 
ment of the superfluid helium film thickness. Our ex¬ 
periment shows a strong dispersive shift due to the elec¬ 
trons that is many times the cavity linewidth. On aver¬ 
age, the coupling per electron in the ensemble is about 
1 MHz, suggesting single electron strong coupling should 
be within reach. Electrons can be held for many hours, 
and their normal mode frequencies and number can be 
controlled by adjusting the trapping potential. The re¬ 
sulting evolution of the dispersive shift agrees excellently 
with our numerical model. 

EXPERIMENTAL SETUP AND DETECTION 
TECHNIQUE 

The electron on helium circuit QED setup consists 
of an integrated electron trap and coplanar waveguide 
(CPW) resonator (see Fig. la). The ground planes of 
the resonator are thicker than the center pin, forming 
a micro-capillary channel which determines and stabi¬ 
lizes the superfluid helium film thickness [26, 27]. The 
electrons are held in the resonator volume by a DC bias 
voltage applied to the center pin as shown in Fig. lc. 

The electrons are confined in both the transverse and 
longitudinal direction of the channel. In the transverse 
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direction, a DC-voltage applied to the center pin cre¬ 
ates a parabolic trapping potential (Fig. lc) that con¬ 
fines the electron ensemble in the channel, colocated with 
the microwave field. In addition to the large electron 
trap formed by the resonator center pin, the devices also 
contain smaller /im-size electron traps positioned near 
the voltage maxima of the fundamental mode for fu¬ 
ture single-electron experiments (Fig. Id). Those smaller 
traps were set to ground potential throughout the exper¬ 
iments discussed in the rest of this paper. The input 
and output of the resonator (Fig. lb) side of the coupler 
are held at 0 V to prevent electrons from leaking out 


the sides. Along the cavity DC bias lead, where the po¬ 
tential may be positive, constrictions shield the electrons 
ensuring that there is a potential barrier for escape. 

The bare cavity resonance frequency is u> o/27r ~ 4.789 
GHz, loaded quality factor Q L ~ 17750 and correspond¬ 
ing decay rate k/2tt ~ 270 kHz in the absence of any 
superfluid or electrons. The Q of the sample is set by the 
couplers, not by the internal Q of the resonator, despite 
the fact that the DC bias lead directly connects the cen¬ 
ter pin to a low impedance. This is possible because the 
connection is made at a voltage node, where radiation is 
minimized [28, 29] 



FIG. 1. Device, circuit schematic and trap geometry, a) Op¬ 
tical and SEM images of a cavity-electron ensemble trap on a 
2 x 7mm superconducting chip. The device is positioned 5.5 
mm above the bottom of a cylindrical superfluid reservoir of 
radius r = 3.175 mm, mounted in a hermetically-sealed cop¬ 
per box at 25 mK in a dilution refrigerator, b) Interdigitated 
gap capacitors with gap width 2 /im at the cavity input, c) 
DC bias electrode connected directly to the center pin of the 
cavity at a node of the standing wave voltage distribution of 
the fundamental mode, d) Sub-pm-size electron trap near 
voltage maximum of the fundamental mode with constriction 
of width 500nm. e) Circuit schematic showing the voltage 
distribution of the fundamental mode (red) and the simpli¬ 
fied measurement and control circuit connected to the center 
pin (pink). The cavity is measured in transmission using a 
low-noise amplifier and the trap potential is tuned through 
a DC source connected to the center pin (pink) through a 
low-pass filter, f) Cross-sectional view of the cavity waveg¬ 
uide gap showing the schematic trap geometry. The ground 
planes (gray) form a micro-channel of height d = 800 nm 
and width wg = 6 /im filled with superfluid 4 He by capillary 
action. A DC voltage on the submerged center pin (pink) 
of width wcp = 2 fi m and thickness t = 80 nm creates a 
parabolic trapping potential for electrons above the surface 
which couple to the RF field in the cavity. 


HELIUM DYNAMICS 

An important prerequisite for trapping electrons on 
helium in a micro-channel geometry is to establish a self- 
stabilized film of superfluid helium of known thickness, 
which can be achieved by capillary action filling of the 
channels from the low-lying bulk reservoir. The helium 
raises the effective dielectric constant of the waveguide, 
lowering the resonator frequency proportionally to the 
thickness, h. 

To measure the cavity response to superfluid helium, 
we monitor the resonance frequency and quality factor 
in transmission while increasing the bulk helium reservoir 
level in small increments (the center pin is held at ground 
potential throughout this measurement). The results of 
such a helium filling experiment are presented in Fig. 2 a. 
Four different regimes can be clearly distinguished in the 
frequency shift. For small amounts of superfluid (regime 
I), an unsaturated van-der-Waals (vdW) film of thickness 
h ~ 30 nm forms as the liquid evenly coats the surface 
of the resonator and the interior of the sample cell, lead¬ 
ing to small frequency shifts of Auio(h)/2n ~ —190 kHz. 
Once the vdW film has saturated, the liquid film shape is 
determined by capillary action with a semi-circular pro¬ 
file z(x) ~ x 2 /2R c (H) in the gap. The capillary radius 
R C {H) = < 7 / pgH is determined by the distance between 
the bulk helium level in the reservoir and the chip surface 
H, where a = 0.378 x 10 -3 N/m is the surface tension 
of liquid helium in vacuum, p = 0.154 x 10 -3 kg/cm 3 
the mass density and g the gravitational acceleration. 
When the radius of curvature becomes on the order of 
the gap width R c ~ wg, the gap starts to fill up by 
capillary action and is filled completely for R c wg 
(regime II). Small increases in shift in the subsequent 
“flat” regime are due to decrease in the curvature of 
the helium profile. Finite element simulations show that 
when the channel is filled (h = 800 nm) the frequency 
shift is —8 MHz, in good agreement with the observed 
data. As H —> 0, the radius of curvature becomes on 
the order of the chip dimensions and eventually starts 
to diverge. We attribute the abrupt jump at 170 mm 3 
to the formation of a thick him which spans the entire 









3 


Relative bulk reservoir level, H (mm) 

5.5 4.8 4.1 3.3 2.6 1.9 1.2 0.4 - 0.3 



IV 



03 

3 

O' 


FIG. 2. Cavity response to superfluid helium, a, Measured 
resonance frequency shifts Auiq/2-k (blue dots, left axis) and 
loaded quality factor Ql (green triangles, right axis) as func¬ 
tions of superfluid volume supplied to the cell (bottom axis) 
and relative bulk helium level H in the reservoir (top axis). 
Each datapoint corresponds to an increase in superfluid vol¬ 
ume of AVrf ~ 2.3 mm 3 and reservoir level AH ~ 70 ;im. b, 
different filling state corresponding to the different regimes 
in a. c, Frequency shift (blue dots, left axis) and quality 
factor (green triangles, right axis) as functions of center pin 
voltage bias V cp at fixed helium level in the capillary action 
regime indicated by arrow in a. Gray dots are frequency shifts 
extracted from single-shot cavity transmission measurements 
with IV = 80 such measurements per voltage bias point V cp . 
The blue data points are averages over the single-shot mea¬ 
surements at each point. 


linearly (regime IV) until the resonator becomes insensi¬ 
tive at thickness h ~ 6 /mi, corresponding to a frequency 
shift of Awo(/i)/27t = —14.145 MHz, again in good agree¬ 
ment with numerical simulations that predict a final fre¬ 
quency shift of —14.1 MHz. All subsequent experiments 
are done at the filling level indicated by the black ar¬ 
row in Fig. 2 a, which corresponds to h ss 647 nm where 
the frequency shift is Au!o(h)/2n ~ —7.4 MHz, and is 
only slightly changed by additional fluid introduced to 
the reservoir. 

The superfluid level in the resonator gap can be mod¬ 
ulated in-situ by sweeping the voltage of the center pin 
V cp . To lowest order, the equilibrium film thickness at 
the center of the gap is determined by the electrome¬ 
chanical force on the film surface and surface tension 
with a quadratic voltage dependence h(V cp ) ~ h(0) + 
(17p/16cr)(eHe — £o)- Fig- 2 b shows measured frequency 
shift and quality factor as functions of center pin voltage 
at a reservoir level of H ~ 4 mm where the gap is partly 
filled by capillary action. The resonance frequency shows 
a parabolic voltage dependence while the quality factor 
remains constant, as expected. The maximum observed 
shift at V cp = +6 V of Awo/27t = —100 kHz corresponds 
to a change in film thickness of Ah — 13 nm at the cen¬ 
ter of the gap. The slight offset of the other-wise sym¬ 
metric response is not well-understood, and it is device 
dependent. The frequency sensitivity to level changes 
in the gap can be estimated from finite element electro¬ 
magnetic simulations to be « 8 kHz/nm, consistent with 
the overall slope and frequency shift. The stability of 
the capillary action film is estimated from consecutive 
single-shot transmission measurements (gray data points 
in Fig. 2 b), with N = 80 frequency measurements per 
voltage bias point. Slow fluctuations of the helium level 
are manifested in the /2n = 16 kHz scatter of res¬ 

onance frequencies (gray points in Fig. 2b) corresponding 
to helium level fluctuations of ShS rlns ^ ss 2.6 nm. 

In summary, the microwave measurement provides a 
high bandwidth way to measure the helium level and its 
fluctuations, down to pm/-\/Hz level sensitivity. Using 
this technique, we establish a lithographically defined, 
stabilized superfluid helium film within the cavity-trap. 
In the regime of the experiment, the channel helium level 
is insensitive to the small differences in the amount of he¬ 
lium put into the sample box. The measured helium level 
fluctuations are relatively small and will be monitored to 
see if they have a significant effect on the electron coher¬ 
ence time through changes in the trapping potential. 


chip (regime III) and is supported by the sample holder. 
As R c -¥ oo, the superfluid film becomes sensitive to me¬ 
chanical vibrations and small fluctuations in the reservoir 
level which manifests itself in a perceived drop in quality 
factor in this regime. Once the reservoir has been com¬ 
pletely filled, the helium fills the region above the chip 


DISPERSIVE MEASUREMENTS OF 
ELECTRONS IN A CAVITY 

Following uncharged superfluid measurements, we pro¬ 
ceed to load electrons into the cavity mode volume 
and detect the trapped electron ensemble in transmis- 
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FIG. 3. Detection of a trapped electron ensemble on super¬ 
fluid helium in a cavity transmission experiment, a, Nor¬ 
malized transmitted power through the cavity as a function 
of center pin trap voltage V cp . b, Normalized transmission 
spectra at V cp = +0.91 V and +2.8 V, showing a shift in 
resonance frequency and a reduction in transmitted power at 
the bias points indicated by the dashed vertical lines in a. 
Solid red lines are fits to Lorentzians. c,d, Resonance fre¬ 
quency (c) and loaded quality factor (d) as functions of trap 
bias in the presence (blue) and absence (red) of an electron 
ensemble. In Figs, a - d, electrons are first loaded into the 
cavity mode volume at an initial bias of V cp = +3V and a 
fixed helium level in the capillary regime with an uncharged 
shift of Acuo/27t = —7.58 MHz and a reservoir level of H ~ 4 
mm. The blue line shows the trap voltage being swept from 
+3 to -IV and back in 4 mV steps, eventually depleting the 
trap region, while the red line shows the same sweep for an 
empty trap. 


sion measurements. Electrons are generated via pulsed 
thermionic emission from a tungsten filament mounted 
in vacuum above the device and attracted towards the 
superfluid surface in the resonator channel by a posi¬ 
tive trap voltage V cp . After waiting for the sample to 
cool, the cavity transmission is monitored while tuning 
the center pin voltage V cp (Fig. 3) starting from +3V. 
The dispersive interaction of the cavity with the trapped 
ensemble leads to a voltage-dependent shift of the cav¬ 
ity resonance towards lower frequencies before reverting 
back at negative trap potentials. We observe maximum 
resonance shifts of up to Ac j > 10k cavity linewidths in 
frequency while Q is somewhat reduced (Fig. 3 b). The 
electron-induced frequency shift reaches a maximum of 
Acu max /27r = —2.47 MHz at V r C p h ' ) = +0.91 V with a 
drop in quality factor and a corresponding increase in 
cavity decay rate of AK max /27r = 122 kHz (blue curves 
in Fig. 3 c and d). These changes in the cavity resonance 
frequency are at least an order of magnitude larger than 
those caused by the electric field induced helium film 
thickness change without electrons. Below the thresh¬ 


old V,-p b ‘ ! , the electron-induced shifts decrease gradually 
as electrons are lost from the trapping region. To ensure 
that the observed cavity response is due to the trapped 
electron ensemble, we performed a control experiment 
where the filament was fired while the center pin was 
biased at -IV. The voltage is then swept in the reverse 
direction (red curves in Fig. 3c and d). The voltage de¬ 
pendent signal is completely absent (see Fig. 2 b and dis¬ 
cussion above). The ensemble-induced cavity response 
has been reproduced in independent experiments using 
five different devices. The maximum observed resonance 
shifts are repeatable and generally vary between 2 — 8 
MHz based on loading conditions. 



FIG. 4. Measured cavity frequency shifts as a function of 
trap voltage bias 1+ P - Colors indicate consecutive cycles of 
the voltage sweep. In each cycle V cp is decreased (arrows) 
until electrons are irreversibly lost from the trap and then 
increased to the initial value of 3.0 V. Gray dashed lines in 
the background show frequency shifts predicted by molecular 
dynamics simulation. For each of the four iso-electron number 
curves the number of electrons in the simulation is depicted 
to the left. A red dotted line indicates modeled electron loss 
with respect to a leak voltage of 530 mV. The iso-electron 
number curves from the experiment terminate within a small 
neighborhood of this loss frontier. 

To further investigate the quantitative form of the cav¬ 
ity shift as a response to the number of electrons N and 
the bias voltage V cp , we designed a protocol to partially 
drain the electrons from the trap. In the following experi¬ 
ment at the end of the first sweep from +3V (bottom half 
of blue curve in Fig. 4), we deliberately sweep down be¬ 
low the V C p h) to introduce electron loss. Now with fewer 
electrons in the trap (top half of blue curve in Fig. 4), the 
cavity resonance shift is smaller in magnitude than be¬ 
fore, but non-hysteretic unless another electron-loss event 
is triggered. For subsequent sweeps, we then set succes¬ 
sively lower stop voltages. The final sweep (orange curve 
in Fig. 4) depletes the mode volume of all electrons as 
evidenced by the vanishing frequency shift on the final 
upward sweep. By carefully tuning the stop voltages and 
number of sweeps we can controllably reduce the number 
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of electrons as desired. 

A quantitative understanding of the electron-cavity in¬ 
teraction requires modeling of the classical many-body 
interactions between electrons as well as the coupling of 
the resulting electron normal modes with the cavity. We 
develop a non-perturbative numerical model which de¬ 
termines the electron ensemble configuration, frequencies 
and coupling to the cavity. First simulated annealing is 
employed to find the minimal energy configuration of the 
electrons. The electrostatic potential used in this step 
is constructed using field profile data derived from the 
sample geometry. After obtaining the electron configu¬ 
rations, we solve the equations of motion of the cavity- 
electron coupled system in a non-perturbative way to 
obtain the cavity frequency shift. For a given number 
of electrons in the trap, this calculation is repeated for 
various I4p to produce the corresponding frequency shift 
curve. Curves for various number of electrons are com¬ 
puted, and no fitting parameters are used in the model 
besides picking the closest iso-electron number curve. 
Fig. 4 shows excellent agreement between the data and 
our computational model. The model is described in de¬ 
tail in the supplementary material. 

Using the measured signal and our model, we can in¬ 
fer the electron mode frequency and coupling strength g. 
The frequency of the strongest coupled mode grows pro¬ 
portionally to ^yVcp and is roughly 25 GHz at V cp = 0.5 
V, which shows we are well within the dispersive regime. 
In this regime the coupling is linearly proportional to 
the distance of the electron to center of the trap, approx¬ 
imately 1.6 MHz//im per electron. In future trap designs, 
additional guard electrodes can improve trap stability at 
low voltage allowing access to the resonant regime. The 
coupling though quite large already, can be enhanced fur¬ 
ther by shrinking the dimensions of the trap. 

There are two types of electron loss observed. The 
first, which determines the number of electrons loaded, 
occurs at higher V cp due to hydrodynamic instability [26]. 
The density immediately after loading at V cp = 3V was 
n ss 2 x 10 9 cm -2 . The second type of electron loss oc¬ 
curs when the trap depth becomes sufficiently shallow, 
such that electrons can leak out of the trap. We model 
this phenomenologically by assuming that electrons are 
lost if the potential difference between ground plane and 
electrons is less than Ui ea k- Using Ui ea k as a single fit 
parameter in the molecular dynamics simulation, we find 
best agreement between simulation and experiment when 
Ueak = 530 mV. The residual population for small elec¬ 
tron numbers at small V cp (orange curve in Fig. 4) is not 
well understood. 

In summary, we have demonstrated the successful trap¬ 
ping and detection of an electron ensemble above the 
surface of superfluid helium in the circuit QED architec¬ 
ture. The measurement technique introduced here could 
extend traditional electrons on helium experiments to 
smaller ensembles and enable observation of the electron 


dynamics. The observation of the large dispersive shift 
and the good agreement with our numerical simulations 
indicate that it should be possible to perform cavity QED 
experiments in a single electron quantum dot. Though 
small, the fluctuations in the helium film thickness are an 
important source of decoherence for the electron motional 
states, and merit further study. Finally, the sensitivity of 
the device to helium thickness changes can be exploited 
for novel cavity optomechanics experiments with super- 
fluid ripplons. 
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Nanoscale Materials, an Office of Science user facility, 
was supported by the U. S. Department of Energy, Of¬ 
fice of Science, Office of Basic Energy Sciences, under 
Contract No. DE-AC02-06CH11357. 
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Supplemental information 

The situation where many interacting electrons in an electrostatic trap couple to a microwave resonator represents 
an interesting physical system. Experimentally it is not possible to directly observe how electrons arrange themselves 
inside the trap. However, such configurations contain important information of the electrostatic and microwave 
properties of the system. Therefore, the goal of this supplement is to develop a set of tools that provides a complete 
quantitative understanding of the classical interaction between few electrons in a trap and a microwave resonator. 
The first step is to start with an electrostatic simulation and solve for the equilibrium position of each individual 
electron in a given trapping potential. 

With the electron configuration at hand, we can then calculate the dispersive shift of the cavity. We take a non- 
perturbative approach to directly compute the cavity frequency from the equation of motion of the electrons-cavity 
coupled system. The numerical results are confirmed by a simple analytic model in the low density regime. In the high 
density regime where our experiment lies, the simulation model agrees well with the experimental data. Using this 
model, we are able to derive a number of relevant quantities, such as the total number of electrons on the resonator, 
and the average coupling per electron. 

The fabrication recipe and experimental setup can be found at the end of the document. 


Electrostatic Simulation 


The first step of our modeling effort starts with the electrostatic simulation of the trapping potential. We simulate 
the trapping potential for a CPW in Maxwell, a finite element simulation tool. Since the geometry does not change in 
the y-direction, we only extract the cross-sectional profile. This numerical potential profile is then fitted with a 14 th 
order polynomial (with only the even power terms). From this we extract the second order coefficient to construct 
an ideal parabolic potential for use in the molecular dynamics simulation. Our analysis shows that for the most part, 
the electrons stay within 500 nm from the center, where this approximation is good. 
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FIG. SI. Potential profile of the trap, showing the original numerical profile derived from a finite element simulation model; a 
polynomial fit,and a parabolic approximation using the second order term of the polynomial fit. 

Numerically, the electrostatic potential used in the molecular dynamics simulation is 

U{x,V c p ) = 0.0733 eU cp ;r 2 , (SI) 

where x is measured in fj ,m from the center of the trap and V cp is the bias voltage applied to the center pin. 
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Molecular Dynamics 


To solve for the equilibrium positions of the electrons in a given electrostatic potential, we use a simulation method 
called molecular dynamics. Molecular dynamics represents a class of deterministic, iterative algorithms that can be 
used to find solutions to many-body problems. For this work, we use the package HOOMD [30]. This general-purpose 
particle simulation toolkit scales on multi-core CPU and GPUs which allows us to quickly anneal up to 4000 electrons. 

In a typical simulation run, we start with a fixed number of electrons n. Due to the large number of pairwise 
interactions that grows quadratically with n, we shrink our 12 mm long resonator from the actual length down to a 
50 /im x 50 n m box with periodic boundary conditions on each side. Simulation and experimental results can then be 
compared by multiplying the number of electrons in the 50/cm x 50 /mr box by L res /Lb ox ~ 243. In the remainder of 
this supplement we refer to the number of electrons on the resonator as N = 243 • n. 

There is no long-range screening in our simulation. However, the Coulomb interaction is cut off at 20 /im to prevent 
electrons from interacting with their own image charges across the periodic boundary condition, causing an explosion 
of pairwise interactions in the system. 

The lowest energy electron configuration is found by annealing the system. Here the temperature of the ensemble is 
gradually decreased until a temperature of below 1 K is reached. Starting from a random initial electron distribution, 
a typical annealing procedure for n = 1000 (N = 2.43 x 10 5 ) in a 50/cm x 50 /im box takes 20 minutes. The 
simulation has two input parameters: N and U cp , the second of which determines the depth of the trapping potential. 
We simulate several N 1 s for the range of U cp used in the experiment. Results from multiple runs with different random 
initial conditions were consistent. 

Fig.S2 shows equilibrium electron configurations for six different pairs of {N, I4p}- The electrons arrange themselves 
in rows along the //-direction, the direction perpendicular to the electrostatic trap. Increasing the trap bias voltage 
leads to an increase in electron density, a reduction in the number of rows, and a decrease of the overall width of the 
ensemble. This is depicted more clearly in Fig. S3 for N = 2.43 x 10 5 . In this figure the color represents the binned 
electron density along x (bin size 6 nm), which was obtained by integrating the electron distribution along y. Gradual 
transitions from 9 rows to 4 rows can be observed as U cp is swept from 0 to 4V. 
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FIG. S2. Examples of electron configurations obtained by molecular dynamics simulations, (a)-(c) Equilibrium electron 
configurations for N = 0.61 x 10 5 electrons on the resonator as the bias voltage is increased from 0.5V, 1.5V to 3.0V, respectively. 
(d)-(f) Same as in (a)-(c) but for N = 2.43 x 10 s . 
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FIG. S3. Integrated electron distribution for a large number of electrons on the resonator: N = 2.43 • 10 5 . For each bias 
voltage the density along the a~-axis was obtained by integrating the electron distribution along y. Gray dashed lines indicate 
transitions that lead to configurations with one fewer row. 


Equation of motion for electrons coupled to a CPW cavity 

With the electron configurations in hand, our goal is to compute the dispersive resonance frequency shift. One 
approach is to calculate all the electron normal modes, and use perturbation theory to find the frequency shift. This 
approach is difficult, as defects in the electron Wigner crystal and the flat potential along the y-direction lead to 
degenerate modes. To tackle this problem, we follow a non-perturbative approach, and calculate the cavity normal 
mode from the equations of motion of the entire electron-cavity coupled system. 


Electron Subsystem 


To calculate the equations of motion for the electrons, we use a Lagrangian formalism. The full Lagrangian C 
consists of an electronic part C e , a cavity part C c and an interaction energy ^coupling- The electronic part is given by 


Ce 



-e^V D c(ri) 
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1 


2 47re 0 


EElT-t: 




(S2) 


where r; = ( Xi , yi) is the coordinate of electron i and VDc( r i) is the electrostatic potential that defines the static trap 
for the electrons. In our case Vbc( r ) = ^trap^ 2 - The trap depth and thus fc tra p are determined by the center pin 
bias voltage. 


Cavity 


Working in the charge basis, the cavity part of the Lagrangian is 


Cc 



Q 2 
2 C’ 


(S3) 


where L and C are the effective inductance and capacitance respectively. 
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Electron-Cavity Coupling 


The coupling between the electrons and the cavity can be written straightforwardly as 

^coupled = C 'y ' VrF (rj, () j. 
i 


(S4) 


Here Vrf^oQ) is the RF potential generated by the resonator. It should be emphasized that generally this potential 
may have a different position dependence than the electrostatic potential VdcM = However, in this work 

the electrostatic potential is determined by the bias voltage on the center pin, whereas Vrf^jQ) depends on the 
RF voltage (or rather charge) on the center pin. Since both potentials originate from the center pin, the functional 
dependence is the same. 

To proceed with the analysis, the charge and position dependent parts are separated, such that 


C 


coupled 


eQ 

~C 


E H RF (ri). 


(S5) 


t/Rpfrj) is a dimensionless function that describes the position dependence of the microwave potential. 


Full Equations of Motion 

The equations of motion are easily obtained by calculating the derivatives with respect to Xi, yi, Q and their time 
derivatives. In general 


which for the cavity results in: 


d dC _ dC 
d t dq dq 


iQ+§ + ^E C/ RF( r *) = ° 

i 


In a similar fashion, the equation of motion in the ^-direction for electron i becomes 

„2 


m e Xi 


<9Vdc f \ , e Q 9 Urf , x 1 e" 

e ~ -( r ‘) + ( r 0 - ol 


dx. 


C dxi 


2 47rei 


E 


Xi — Xj 


= 0. 


(S6) 


(S7) 


(S8) 


Eqs. (S7) and (S8) are the exact equations of motion but are also highly nonlinear. To this end the system is 
linearized around the equilibrium point (:r^ eq , Vi,eqi Qeq = 0). This means that the equations of motion will contain 
terms involving r ieq . This underlines that obtaining the equilibrium electron configuration via molecular dynamics 
simulations is of critical importance. 

At the equilibrium point the potential energy 


W= ^ + ^E C/ RF( r 0 + eE^c(ri) + JtE-EE^ 


2 47re 0 ^ ' |r* - r^. 


of the system is minimized, such that 
dU 


dx~ Vi,eqi Qeq — 0) ^ (l\,eq) 

^L( x . v . Q — o) — e—^2-(r- ). 

Qy v^^eqj yi,e q? ^eq — v j — 0 Qy 


f ^ %i,e q *^f,eq q 

2 47reo |r ljeq -r iieq | 3 


1 e 2 


E 


Vi,e q Uj,e q 


2 47re 0 |r,;. fin - r,, 


= 0 


L 2 ,eq A j,eq| 


~QQ(%i,e qj Ui,e qi Qeq = 0) = — ^ ' C7 R p(r,; jeq ) = 0 


(S9) 
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(511) 
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Finally, Eqs. (S7) and (S 8 ) can be linearized. First consider Eq.(S7). If 5Q is a small deviation from Q eq = 0 and 
Sxi is a small deviation from Xi, eq this equation may be written as 

L6Q + ~C + C E 

i 

and by using Eq.(S12) we arrive at the equation of motion for the cavity: 

L6Q + ~c + c E 


9Urf 

dx i 


(U.eq)*^/ 


9Urf 

dyi 


(r it eq)6yi 


= 0 ,. 


(S14) 


TT , \ , 9Urf , . 9Urf , 

CRF (U,eq) T (U,eq)oX/ -f- „ (U,eq)cZ /2 


dxi 


dyi 


= 0 , 


(S13) 








FIG. S4. Schematic representation of two nearby interacting electrons located at (xi,yt) and ( Xj,yj ). The distance between 
the two electrons is given by | — ry | = rij. The angle between the line that connects the two charges and the x-axis is dij. 

Here the trap confines electrons along the x-axis. The electron spacing along the {/-direction is ' and finally the row spacing 
is d. Note that this figure is not to scale. 


Now for the electrons we linearize Eq.(S 8 ). To do this, it is useful to note that (see Fig.S4) 
Xi — Xj rij cos dij + (Sxi — 8xj ) 


ij 

cos 9 


'«* »«(i + (i + iSS/T 


-| 3/2 


ij f | q /rift \\8xi Sxj 3 . . 5yi Syj 

tE ~ b ( 1 + 3cos(2%))- -3 3 - - - sm( 2 %)- 

* j * ' ij ^ ' *i 


0(S 2 ). 


Here = |r/ ieq — ry eq | and Oij is the angle between r,; eq — ry e q and the x-axis. 
The linearized version of the equation of motion then becomes 


m e Sxi + e 
1 e 2 


, , d 2 V DC , , d 2 V DC , ^ 

dxi (ri ’ eq)+ dx 2 {ri ’ eq )‘ 5Xi+ 5x i % (ri ’ eq)% 


e5Q dUR f 

~c feT 


( r i,eq) 


2 47re 0 ‘— d . r 


Using Eq.(SlO) this can be simplified to 


_ 1/1 3 

y, ( Tij cos Oij - - (1 + 3cos(2%)) ( Sxi - Sxj ) - - sin(2 %)(% - Syj) 


= 0 


m e Sxi + e^uKeq^x* + e^^-{v l ^)5y l + ^^^(r heq )SQ 


1 e 2 


4 47T£o 


E 


(1 + 3cos(2%)) SXi ~ +3sin(2%) % ~^ 

r b ' b 


= 0 . 


(S15) 


(S16) 


(S17) 
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Finally, for the sake of transparency let 


;± 1 e 2 l±3cos(2 Oij) A 1 1 e 2 3sin(2 0^) 

/Cl y — " n ana to n — ~ 5 

- 3 4 47 T 60 rf. 


ij 4 47T60 
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such that Eq.(S17) simplifies to the linearized equation of motion for the electrons: 


(SIS) 


Sx i + A?^( r *.eq )*Q + I et 7)3 2 ( r *.eq) + E k % 1 Sx 


d 2 Vi 


C dx.j 


dx 2 




E k tj^ x 0 + ( e Qx Qy. ( r *> e q) ^ I ^/* E kjfyj — 


(S19) 


Above equations show that a change in position (either in x or y) of electron j translates into an x displacement of 
electron i. In the y-direction the equation of motion follows in a similar way: 


m e 8yi + ^^^-( rj ,eq )5Q + I e U y ° c (r,, eq ) + ^ I Sy, 


d 2 V T 


C dyt 


dy 2 




E k ij^yi + ( e dx dy ^ r * ,eq ^ I E kj^ x j ~ 

V 1 1 j 


(S20) 


The system of equations (S14), (S19) and (S20) can be written in matrix form as follows: 


/SQ\ 


(SQ\ 
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V i ) 


\ i / 
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(S21) 
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With the kinetic matrix and mass matrix, the eigenvalue problem 

=Ui\r)i) 
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(S24) 

can be solved to obtain the normal modes \rji) and corresponding frequencies w*. The eigenfrequency corresponding 
to the mode with the highest cavity participation tells us the dispersive cavity shift. 


Constraining the Motion Along the ir-axis 


In the previous section the most general equations of motion for a two dimensional confined electron gas were 
derived. In this case, the static trap is flat in the y-direction: 

Vbc(r») = ^fctrap£?, (S25) 

and Us. f has the same r-dependence as Vdc : 

UnF^i) = Pxl- (S26) 


Since both Vbc and 1 /rf do not depend on yi the equations of motion simplify drastically. Additionally, since the 
trapping potential is flat along the y-direction and due to the presence of defects - either intrinsic or due to imperfect 
annealing - one can get modes with imaginary or zero frequency. To reduce the impact of these modes we “freeze” 
the modes in the y-direction, i.e. we set 5yi = 0 in what follows to simplify the equations of motion even further. 

After these simplifications the matrices that govern the equations of motion read 


and 


(L 0 0 ...\ 

0 m e 0 ... 

0 0 m e ... 


V ; ; ■■■/ 
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2 GpX\ >e q 
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2e/3x 2l eq 
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2e/3iri,e q 

2e0X2,eq 

2e/3x 3i a q 
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£^trap "f" j 

— fcl2 

-ki3 

— ^21 

^trap ^2 j 

^23 

-k 3 1 

—k 3 2 

efctrap + ksj 


V 




(S27) 


(S28) 


Analytical Model of Simple Electron Configurations 


With only two rows of electrons inside the channel, it is possible to find the lowest energy configuration analytically. 
This serves as a good comparison for our numerically simulated electron configurations. Let the transverse direction 
of the trap be x, then the trap is parabolic along the x direction and flat along y. If we assume no electron loss 
and only two rows, as we vary the trap bias, only the transverse configuration of the electrons changes. The electron 
density along the y axis (A) remains constant. 

In the simple geometry depicted in Fig.S4, the force electron i exerts on its neighbor j along x is 


Fn 


x , single neighbor 


1 e 2 d 
47 re 0 r?. r i:j ' 


(S29) 


Since electron j has another neighbor to the right, the resulting force is close to 2 F x , ee . If we take into account of the 
next-nearest neighbor, the total force one electron experiences is 


4/;.total — 'yFx, single neighbor; 


(S30) 
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where 7 ss 2.08. Depending on the electron density, 7 grows from 2.08 in the sparse limit to around 2.38 in the dense 
limit where the ensemble transits into 3 rows. 

Now to find the equilibrium configuration, we have 


-F s 


x,trap 


= 7 Fx, 


single neighbor’ 


(S31) 


Note that the x -location of the electrons is measured from the center line, therefore x = d/2. We can then rewrite 
Eq.(S31) as 


ek. 


trap 


7e 


2 47reo (l/A ) 2 + d 2 ri 


(S32) 


Here fc trap is the curvature of the trap, which grows linearly with respect to the trapping bias voltage I4 p . For any 
bias voltage, we can now solve for d in the equilibrium configuration 


d z = 


7 




2/3 


Hrap 


27T£o 


- (1/A ) 2 


(S33) 


In Fig. S5 we plot the width of the ensemble for different densities A. 



FIG. S5. 2-row ensemble width as a function of trap bias voltage V cp at various electron densities. An electron density of 
A =0.8 (2.0) /jm" 1 corresponds to N = 9.7 x 10 3 (N = 2.4 x 10 4 ) electrons on the resonator. 

For simple configurations that have only 1 ~ 2 rows, we derived the analytical solution of the electrostatic equilib¬ 
rium configuration. However, as the electron density and the number of rows increase, the total number of interactive 
terms that need to be taken into account quickly grows to become unmanageable. More importantly, because the 
trapping potential is more shallow on the side than in the middle, adjacent rows of the electron ensemble are not 
commensurate. This means when there are more than 2 rows, the electrons cannot form a perfect crystal in a chan¬ 
nel with parabolic transverse profile and our analytical solution breaks down. Therefore, for more than 2 rows we 
calculate the equilibrium position of the electrons numerically using our molecular dynamics formalism. 


Modeling Electron-Induced Cavity frequency Shift 

Using the equation of motion of the constrained electron-cavity coupled system, we can now calculate the electron- 
induced cavity frequency shift at various bias voltages V cp with different number of electrons N. First, we calculate 
the normal mode frequency of a simple 2 -row electron ensemble, using the analytical solution to the equilibrium 
configuration (Eq.(S33)). Then we compare this result with the normal mode frequency calculated from electron 
configurations obtained using HOOMD. The results are depicted in Fig.S 6 . 
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FIG. S6. Comparison between the analytical model and the molecular dynamics simulation for two different numbers of trapped 
electrons: N = 0.41 • 10 5 (left column) and N = 0.61 ■ 10° (right column). Top: binned electron density (bin size = 6 nm) as a 
function of x via HOOMD simulation (color in the background) and calculated width of a two row electron configuration obtained 
from the analytic model (blue). Note that there is no free parameter. Bottom: Cavity frequency shift solved numerically (with 
HOOMD, orange dots) and analytically (blue). In the analytical solution, a simplified version of the equations of motion was 
used. 


In Fig.S 6 a and S 6 b., we compare the ^-locations of the electrons in the ensemble from HOOMD with those from 
the analytic model. In the analytical model, we only take into account nearest neighbor interactions, such that the 
geometric parameter 7 in Eq.(S33) is equal to 2.08. Clearly, this results in good agreement in the region where 
V cp > 1.0V. For lower bias voltages it is energetically more favorable to form three rows instead. 

Next we compare the cavity frequency shift for both cases. The cavity frequency shift is obtained by solving 
the equation of motion, Eq.(S24). For the analytical model we exploit symmetry to simplify the kinetic matrix 1C. 
Therefore, the problem reduces to diagonalizing a 2 x 2 matrix A4 _ 1 A, where 


/ u 2 0 L depu%Ly/N\ 
\dej3ulL\fN e/c tr ap + 2fci 2 / 


(S34) 


Here wo = y/l/LC is the resonance frequency of the resonator and f3 = 0.0733 m 2 from Eq.(Sl). 

The results are shown in Figures S 6 c and S 6 d. Again, both plots show good agreement between the analytical 
model and molecular dynamics simulation. The largest deviation occurs when the bias voltage is low, where the 
electron configuration transits from 3 to 2 rows. As expected, the 2 -row analytic model no longer describes the 
correct geometry in this regime. 
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Comparison with experimental data at higher densities 


The electron ensembles encountered in our experiment are much more complicated than the simple cases mentioned 
above. However, using the molecular dynamics simulation package HOOMD [31], we are able to anneal systems far 
beyond the simple 2-row case, and solve for the cavity frequency shift using these configurations. 



FIG. S7. Comparison of the measured cavity shift (gray) to the simulated cavity shift (color). The number of electrons in the 
simulation was increased from N = 0.49 x 10 5 (blue curve) to N = 2.43 x 10 5 . For the intermediate curves the number of 
electrons N (in units of 10 5 ) is depicted next to each curve. 

The key observations in our experiment are the smooth change of the cavity frequency shift during the bias voltage 
sweep, and the sharp jumps corresponding to irreversible electron loss from the trap. In our numerical model, we are 
able to replicate both of these phenomena by calculating the cavity frequency shift for each point in our V cp sweep, 
and for various number of electrons N in the trap. For each curve in Fig. S7, we do not rescale the x or y axis nor is 
there any offset applied. The only free parameter is N. By comparing the simulation with the experimental data, we 
are able to estimate how many electrons are present for each trace in the experiment shown in the main text. The four 
curves that match the experimental data have N = 0.61 x 10 5 , N = 1.09 x 10 5 , N = 1.70 x 10 5 and N = 2.43 x 10 5 . 

In the experiment, an important observation is the loss of electrons from the trap. With the equilibrium positions, 
such a process can be simulated. Let us assume there is a leak in the trap with a threshold voltage Ideak that is 
independent of the trap bias. If the screened potential (Fig. S8) that an electron feels is higher than this threshold 
voltage, it is energetically more favorable for this electron to leak out of the trap. Hence the configuration is unstable. 
To determine the threshold voltage, we compute the screened potential for all electrons in each configuration. Then 
we find the threshold voltage that produces the best fit to the experimental data. 

The actual value of the leak voltage determines the starting point of each of the curves in Fig. S7. A leak voltage 
of 530 mV gives the best fit to the data, shown in Fig. S9. 


Normal modes of the electron configuration and coupling strength 

The same normal mode solution to the electron-cavity coupled system also gives access to the relevant modes in the 
electron subsystem. By picking the ten most strongly coupled electron modes, we can estimate the electron normal 
mode frequency during the voltage sweep. 

In Fig.SlOa we plot the normal mode frequency of the trapped electrons. For most I4 p the electron mode frequency 
is tens of GHz, indicating that we are working in the dispersive limit. In this limit the cavity shift is given by Ngf ms / A, 
where g rms is the rms electron-cavity coupling and A is the frequency difference between the electron normal mode 
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FIG. S8. Comparision between the trap potential (light blue) and the screened potential when electrons are present (red). 
A green rectangle in the middle shows where the electrons appear. The rugged shape of the screened potential is due to the 
discreteness of the electron configuration. 



FIG. S9. Loss frontier calculated using V[ ea k = 530 mV. The model fits well with the hysteresis data until the electron ensemble 
is reduced down to mostly a single row in the trap. 


and the cavity resonance. A closer look at the electron mode evolution for TV = 0.61 x 10 5 reveals an interesting 
feature. For this curve, the electron density is low enough to support only 1 or 2 rows of electrons. As we increase the 
bias voltage, the electron normal mode initially grows as yj I4 P , but rolls off as the angle 9 between each electron and 
its nearest neighbor approaches 7r/2 (see Fig.S4). At 9 = 7r/2 the electrons form a single row such that the coupling 
to the cavity vanishes. The voltage at which this happens is determined by the electron density A and can be modeled 
precisely with our analytical model. 

For higher TV, where there are more than 2 electron rows, the electron normal modes evolve mostly as yJV^p- Minor 
row-reconfigurations occur, which result in small jumps in the normal mode frequency. 

Using the fact that we are in the dispersive limit, we can now look at the rms coupling per electron at different 
V cp and different TV. This is depicted in Fig.SlOb. In general, the coupling per electron decreases as the bias voltage 
increases and higher TV leads to a higher overall coupling. Upon further analysis, we found that this trend is mostly 
due to the linear relationship between the location of an individual electron and its coupling to the cavity. In our 
normal mode solution, for an electron located at a equilibrium position x^eq, the coupling term between the electron 
and a single photon excitation in the cavity is 


2 e/3a: 


i,eq 


c 


If we plot the rms coupling per electron as a function of the overall ensemble width, the relationship is roughly linear 
for constant TV, while ensembles with more rows have a smaller rms coupling. 


Experimental Setup 

The measurements are done dispersively by measuring the shift in the cavity resonance frequency. 
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FIG. S10. (a): Electron normal mode frequency for N = 0.61 x 10 5 (orange), N = 1.09 x 10 5 (red), N = 1.70 x 10 5 (green) and 
N = 2.43 x 10 5 (blue), (b): rms coupling per electron for different values of N. Same color coding as in (a), (c): rms coupling 
per electron plotted as a function of the electron configuration width. Same color coding as in (a). Small gaps in the solution 
appear for V cp that correspond to row transitions. For these points the unstable equilibrium configuration results in solutions 
where a small fraction of the electrons has a disproportionate displacement. 


An Agilent network analyzer is connected to the input and output ports of the cavity with attenuation along the 
input line for photon thermalization and amplification on the output. Low-pass filters filled with the Eccosorb epoxy 
are used to filter out thermal photons that would otherwise enter the sample box and affect the cavity Q. 




FIG. Sll. Left: microwave measurement setup, right: gas handling system. 

To feed helium into the hermetic sample box, particular consideration was taken to thermalize the room temperature 
helium gas with each stage of the fridge. Stainless steel capillary tube was wrapped around and soldered to copper 
cylinders at each stage of the fridge. Above 4K, 0.085” ID (1/8” OD) tubes were used to prevent plugs, while below 
4K the capillary tube had an ID of 0.022” (1/16” OD). 

During the experiment, the transmission spectrum of the cavity is monitored using a network analyzer (NWA). The 
electrons are loaded by pulsing a small filament briefly with a negative bias voltage. As the sample cools down after 
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the loading, we record the cavity resonance frequency and the quality factor. 


Device Fabrication 

The coplanar waveguide resonator chip is fabricated via a two-step e-beam lithography process. The ground plane 
is made of 800 nm thick niobium, whereas the center pin of the waveguide is made of 80 nm aluminum. All patterning 
is done on a JOEL 100 keV e-beam writer. 

Due to the thickness of the niobium layer (800 nm), an aluminum hard mask is required for the etching. Re¬ 
depositioning of aluminum (micro-masking) occurred during the RIE process, which is mitigated by over-etching at 
the end. 

Below we list all the fabrication steps: 

1. First Layer 

(a) Nb deposition and wafer preparation 

(b) Coat 2-inch sapphire wafer with 800 nm Nb 

(c) Coat wafer with aluminum as the dry etch stopper layer 

2. E-beam lithography 

(a) Spin coat ZEP 520 at 3000 rpm for 45 s 

(b) Bake resist at 150C for 180 s 

(c) Pattern with a JEOL JBX9300FS tool; area dose 1600 C/cm 2 . 

(d) Develop ZEP in xylene at 0°C, rinse with DI water. 

(e) Oxygen plasma clean 10 min 

3. Dry etch (reactive ion etching) 

(a) Prepare the chamber by running the A1 etch recipe for 10 min. 

(b) To pattern the aluminum hard mask use BCl3(3.0 seem) and Cl2(24.0 seem) at 20°C, for 1:45 min. 

(c) Prepare the chamber by running the Nb etch recipe for 10 min. 

(d) Then with SFg(25 seem) + Ar (5 seem) at 20°C, for 15 min until etched through. 

4. Second layer: 

(a) E-beam lithography 

i. Spin coat ZEP 520 at 4000 rpm for 45 s. 

ii. Pattern with JOEL tool at area dose of 1600 C/cm 2 . 

iii. Develop ZEP at 0°C, rinse with IPA and DI water. 

iv. Oxygen plasma clean 10 min 

(b) Aluminum lift-off 

i. Evaporate 80 nm of aluminum 

ii. Dip wafer in n-methylpyrrolidone (NMP) heated to 60°C, 

iii. As the resist dissolves away, gently blow away the aluminum and dip into acetone heated to 60°C 

iv. take out and blow dry. 


* David.Schuster@uchicago.edu 

[1] Lyon, S. A. Spin-based quantum computing using electrons on liquid helium. Physical Review A 74, 1-6 (2006). URL 
http://link.aps.org/doi/10.1103/PhysRevA.74.052338http://dx.doi.org/10.1103/PhysRevA.74.052338. 



20 


[2] Platzman, P. M. & Dykman, M. I. Quantum Computing with Electrons Floating on Liquid Helium. Science 284, 1967-1969 
(1999). URL http: //www.ncbi .nlm.nih.gov/pubmed/10373109http: //www. sciencemag.org/content/284/5422/1967. 

[3] Dykman, M. I., Platzman, P. M. & Seddighrad, P. Qubits with electrons on liquid helium. Physical Review B 67, 155402 
(2003). URL http: //link. aps . org/doi/10.1103/PhysRevB . 67.155402. 

[4] Schuster, D. I., Fragner, A., Dykman, M. I., Lyon, S. A. & Schoelkopf, R. J. Proposal for Manipulating and Detecting 
Spin and Orbital States of Trapped Electrons on Helium Using Cavity Quantum Electrodynamics. Physical Review Letters 
105, 1-4 (2010). URL http://link. aps . org/doi/10.1103/PhysRevLett. 105.040503. 

[5] Monarkha, Y. & Kono, K. Two-Dimensional Coulomb Liquids and Solids (Springer, Berlin, 2004). URL http://www. 
springer.com/us/book/9783540207542. 

[6] Shirahama, K. & Kono, K. Dynamical Transition in the Wigner Solid on a Liquid Helium Surface. Physical Review Letters 
74, 781-784 (1995). URL http://link.aps.org/doi/10.1103/PhysRevLett.74.781. 

[7] Shirahama, K., Ito, S., Suto, H. & Kono, K. Surface study of liquid3He using surface state electrons. Journal of Low 
Temperature Physics 101, 439-444 (1995). URL http://link.springer.com/10.1007/BF00753334. 

[8] Williams, R., Crandall, R. & Willis, A. Surface states of electrons on liquid helium. Physical Review Letters (1971). URL 
http://link.aps.org/doi/10.1103/PhysRevLett.26.7. 

[9] Grimes, C. C. & Adams, G. Evidence for a Liquid-to-Crystal Phase Transition in a Classical, Two-Dimensional Sheet of 
Electrons. Physical Review Letters 42, 795-798 (1979). URL http://link.aps.org/doi/10.1103/PhysRevLett.42.795. 

[10] Andrei. Two-Dimensional Electron Systems, vol. 19 (1997). URL http://link.springer.com/10.1007/ 

978-94-015-1286-2. 

[11] Ikegami, H., Tsutsumi, Y. & Kono, K. Chiral Symmetry Breaking in Superfluid 3He-A. Science 341, 59-62 (2013). URL 
http://www.sciencemag.org/cgi/doi/10.1126/science.1236509. 

[12] Kono, K. Electrons on the Surface of Superfluid 3He. Journal of Low Temperature Physics 158, 288-300 (2010). URL 
http://link.springer.com/10.1007/sl0909-009-9962-3. 

[13] Chepelianskii, A. D., Watanabe, M., Nasyedkin, K., Kono, K. & Konstantinov, D. An incompressible state of a 
photo-excited electron gas. Nature Communications 6, 7210 (2015). URL http://www.nature.com/doifinder/10.1038/ 
ncomms8210. 

[14] Papageorgiou, G. et al. Counting individual trapped electrons on liquid helium. Applied Physics Letters 86, 1-3 (2005). 

[15] Rousseau, E. et al. Addition spectra of Wigner islands of electrons on superfluid helium. Physical Review B 79, 045406 
(2009). URL http: //link. aps . org/doi/10.1103/PhysRevB . 79.045406. 

[16] Bradbury, F. R. et al. Efficient Clocked Electron Transfer on Superfluid Helium. Physical Review Letters 107, 266803 
(2011). URL http: //link. aps . org/doi/10.1103/PhysRevLett. 107.266803. 1107.4040. 

[17] Takita, M. & Lyon, S. a. Isolating electrons on superfluid helium. Journal of Physics: Conference Series 568, 052034 (2014). 
URL http://stacks.iop.org/1742-6596/568/i=5/a=052034?key=crossref.2fec45ff72d23fff62247d21b092b69b. 

[18] Gor’kov, L. & Chernikova, D. Concerning the structure of a charged surface of liquid helium. Soviet Journal of Experimental 
. . . (1973). URL http: //adsabs .harvard, edu/abs/ 1973JETPL. . 18. . . 68G. 

[19] Edel’man, V. S. Levitated electrons. Soviet Physics Uspekhi 23, 227-244 (1980). URL http://stacks.iop.org/ 
0038-5670/23/i=4/a=R01?key=crossref.248ala88e99c7bc7faeb6c8f819ebe25. 

[20] Petta, J. R. Coherent Manipulation of Coupled Electron Spins in Semiconductor Quantum Dots. Science 309, 2180- 
2184 (2005). URL http: //www.ncbi .nlm.nih.gov/pubmed/16141370http: //www. sciencemag. org/cgi/doi/10.1126/ 
science.1116955. 

[21] Wallraff, A. et al. Strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics. 
Nature 431, 162-167 (2004). URL http://dx.doi.org/10.1038/nature02851. 

[22] Blais, A., Huang, R.-S., Wallraff, A., Girvin, S. M. & Schoelkopf, R. J. Cavity quantum electrodynamics for superconducting 
electrical circuits: An architecture for quantum computation. Phys. Rev. A 69, 062320 (2004). URL http://link.aps. 
org/doi/10.1103/PhysRevA.69.062320. 

[23] Sommer, W. & Tanner, D. Mobility of Electrons on the Surface of Liquid '{4}He. Physical Review Letters 27, 1345-1349 
(1971). URL http: //link. aps . org/doi/10.1103/PhysRevLett. 27.1345. 

[24] Schoelkopf, R. J. & Girvin, S. M. Wiring up quantum systems. Nature 451, 664-669 (2008). URL http://www.ncbi. 
nlm.nih.gov/pubmed/18256662http://dx.doi.org/10.1038/451664a. 

[25] Devoret, M. H. & Schoelkopf, R. J. Superconducting circuits for quantum information: An out¬ 
look. Science 339, 1169-1174 (2013). URL http://www.sciencemag.org/content/339/6124/1169.abstract. 

http://www.sciencemag.org/content/339/6124/1169.full.pdf. 

[26] Marty, D. Stability of two-dimensional electrons on a fractionated helium surface. Journal of Physics C: 
Solid State Physics 19, 6097-6104 (1986). URL http://stacks.iop.org/0022-3719/19/i=30/a=019?key=crossref. 
C8fb6fb583aac58a288a9a74e9c3faf c. 

[27] Rees, D. et al. Point-Contact Transport Properties of Strongly Correlated Electrons on Liquid Helium. Physical Review 
Letters 106, 1-4 (2011). URL http://link.aps.org/doi/10.1103/PhysRevLett.106.026803. 

[28] Armour, a. D., Blencowe, M. P., Brahimi, E. & Rimberg, a. J. Universal Quantum Fluctuations of a Cavity Mode 
Driven by a Josephson Junction. Physical Review Letters 111, 247001 (2013). URL http://link.aps.org/doi/10.1103/ 
PhysRevLett. Ill. 247001. 1307.2498. 

[29] Petersson, K. D. et al. Circuit quantum electrodynamics with a spin qubit. Nature 490, 380-383 (2012). URL http: 
//dx.doi.org/10.1038/nature11559. 

[30] Anderson, J. a., Lorenz, C. D. & Travesset, A. General purpose molecular dynamics simulations fully implemented on 
graphics processing units. Journal of Computational Physics 227, 5342-5359 (2008). URL http://linkinghub.elsevier. 


21 


com/retrieve/pii/S0021999108000818. 

[31] HOOMD. URL https ://codeblue . umich. edu/hoomd-blue/. 


